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We calculate the optical properties of carbon nanotubes in an external static electric field directed 
along the tube axis. We predict strong Franz-Keldysh oscillations in the first and second band-to- 
band absorption peaks, quadratic Stark effect of the first two excitons, and the field dependence 
of the bound exciton ionization rate for a wide range of tube chiralities. We find that the phonon 
assisted mechanism dominates the dissociation rate in electro-optical devices due to the hot optical 
phonons. We predict a quadratic dependence of the Sommerfeld factor on the electric field and its 
increase up to 2000% at the critical field of the full exciton dissociation. 

Semiconducting carbon nanotubes are direct bandgap materials which have attracted much 
attention recently for nanophotonic applications 1]. Almost fifty years ago Franz and Keldysh 
predicted that a static electric field would modify the linear optical properties of the 3D semicon- 
ductors near their absorption edge P, 0] . They showed that the absorption coefficient decays ex- 
ponentially for photons below the bandgap and shows oscillations for energies above the bandgap. 
The interest in electroabsorption was revived about three decades latter after discovery of the 
Quantum-Confined Stark Effect in 2D quantum well structures P, [||. Large Stark shifts were 
observed in fields directed perpendicular to the 2D planes. In ID carbon nanotubes, excitons were 
predicted to have large binding energies (| and to dominate the absorption spectra [3, 0] , a fact 
which was verified experimentally by two-photon spectroscopy P, [l^ and from the observation of 
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the phonon sidebands in photoconductivity spectra 
nanotubes have interesting scaling properties [8] and they can be as small as those in 2D struc- 
tures and as large as 30% of their bandgap depending on both the nanotube structure and the 
environment. 

It has been long recognized that excitonic effects enhance significantly the electroabsorption sig- 



nal Q. 



An electric field leads to several modifications of the absorption spectrum: 1) modulation 
of the absorption coefficient; 2) growth of the band-to-band absorption spectral weight; 3) shift of 
the absorption peak, known as the Stark effect; and 4) dissociation of the bound exciton. In bulk 
3D semiconductors the binding energy is small and most of the theoretical and experimental focus 
has been on the field induced absorption in the region below the bandgap and on the quantum 
confined Stark effect in 2D structures |l3| |. In carbon nanotubes, the binding energy is large and 
the oscillator strength of the higher lying Rydberg states is infinitesimal, so that relatively large 
changes in the absorption at the first excitonic peak and the first band-to-band absorption are 
expected. 

Here we explore the field induced changes in the absorption spectra of nanotubes for a large 
range of photon energies as a function of tube chirality and dielectric environment. In opto- 
electronic carbon nanotube devices, excitons can be produced by electrom-hole recombination or 
by impact excitation by hot carriers Q, Q, [l(J. Therefore, we also address here the problem of 
the bound exciton dynamics produced in the impact excitation process in an electric field and in 
the presence of the hot optical phonons. We use a Bethe-Salpeter equation solution for an exciton 
in the static electric field directed along the tube axis. Similar to Ref. [8( we employ periodic 
boundary conditions for the unit cell of length L (typically L w 450 nm) where the electron- 
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hole interaction is computed by including both the direct and exchange Coulomb terms 
nanotube diameter dependence of the exciton binding energies obtained by two-photon fluorescence 
excitation spectroscopy [l8] agrees very well with our model calculations Q if we choose e = 3.3, 
a value that we use in the rest of the paper, unless stated otherwise. The external electric field 
potential is added to the total Hamiltonian in the Bethe-Salpeter equation. The potential for the 
relative motion of the exciton has been smoothed out to avoid numerical instabilities: 
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U F (r e - r h ) = eF (z e - z h ) tanh f k I \ - ^IhL 11 , ;|. , 



such that the potential is zero at the unit-cell boundaries at \z e — Zh\ = ±£/2, while at small 
\z e — Zh\ <C L/2 the hyperbolic tangent equals to unity if we choose the smoothing parameter fco to 
be > 30. The absorption spectra are calculated as in Ref.[8[ from the q = excitonic wavef unctions. 

The absorption spectra at different fields are shown in Fig.QJu At zero field, there is no absorption 
in the energy range between the first En exciton at 0.52 eV and the onset of the first band-to-band 
(free electron-hole continuum) absorption An at 0.71 eV. The higher energy Rydberg states have 



0.4 0.5 0.6 0.7 0.8 0.9 

E (eV) 



FIG. 1: Absorption spectra in (16, 8) tube d = 1.7 nm in applied electric fields F in V//im, F — 0.0 
(black), 2.0 (blue), 6.0 (green), 10.0 (red) with excitons (a) and for the band-to-band absorption (b) with 
no excitonic effects. The GW corrections are modeled by the scissor approximation with a scissor value 
equal to the first exciton binding energy of 194 meV, such that the onsets of zero field spectra in (a) and 
(b) coincide. 

an infinitesimal spectral weight in ID, except for the highest energy Rydberg state at 0.70 eV, 
whose energy is in resonance with the onset of the first band-to-band absorption. The interaction 
between this excitonic state and the band-to-band continuum leads to a spectral weight transfer 
to the Rydberg exciton, which appears as a peak in the absorption spectra just below the onset of 
the first band-to-band absorption in Fig. [l}i. 

At a finite electric field, the band-to-band absorption is modulated and grows in intensity, as 
shown in Fig. QJi. Similar to the free electron-hole absorption, in the absence of excitonic effects, 
shown in Fig.[T]D, the band-to-band absorption is modulated with a period proportional to the field 
strength. However, the modifications of the band-to-band absorption in the forbidden region (i. e. 
below An) in the absence of excitonic effects, and correspondingly the absorption below An and 
above the excitonic En peak in the presence of excitons, are very different. In the former case, 
the free particle absorption decays exponentially in the bandgap region as shown in Fig. [TJd, while 
in the excitonic picture an additional peak develops in the forbidden region, as shown in Fig. [Ij,. 
As the field is increased, this peak grows in intensity, moves deeper in the forbidden region, and 
eventually merges with the En exciton absorption peak. We define as the exciton dissociation 
field F c the field at which the first Franz-Keldysh peak coincides with the En exciton energy. 
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At the same time, the "transparent" region between the first two oscillatory peaks broadens and 
the first absorption minimum moves towards higher energies above the zero-field bandgap edge 
An. Above the second E22 exciton (not shown in Fig. QJi), the second band-to-band absorption 
above the second bandgap A22 also shows Franz-Keldysh modulations. However, the absorption 
coefficient does not go to zero at energies below A22, because the first band contributes to the 
total absorption in that region. 
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FIG. 2: The relative band-to-band absorption spectral weight as a function of field in (14, 0) tube with 
d = LI nm (green circles), (16, 0) d — 1.3 nm (blue circles), (f6, 8) d — 1.7 nm (cyan circles), (25, 0) 
d = 2.0 nm (magenta circles), and (24, 8) d — 2.3 nm (red circles). The solid curves are the best fits to 
Eq. © with A = 0.009 and (do, K 3 ) from table E 

The spectral weight growth of the band-to-band absorption with the electric field is due to the 
bound exciton wavefunction mixing with the first band-to-band continuum, which leads to spectral 
weight transfer from the excitonic peak to the latter. As a result, the overall oscillator strength of 
the first band-to-band absorption grows as much as 400% before the excitonic peak starts to overlap 
with the tail of the first Franz-Keldysh oscillatory peak. This corresponds to the field strength of 
about 1/3 of the critical field F c required for the full exciton dissociation. Upon exciton dissociation 
the oscillator strength is being transferred to the band-to-band absorption, which is increased from 
its zero-field value by 10-20 times. The field dependence of the Sommerfeld factor, defined here as 
the fraction of the band-to-band oscillator strength, is shown in Fig. [2] and it can be well fitted by 



5 



(n, m) 


d 


E b 


F 


do 


K s 


Kb 


a 


(13, 0) 


1.03 


281 


169 


0.026 


9.1 


4.2 


6.3 


(14, 0) 


1.11 


288 


198 


0.029 


7.0 


3.0 


9.9 


(12, 4) 


1.15 


277 


179 


0.029 


7.2 


3.2 


7.5 


(16, 0) 


1.27 


236 


119 


0.026 


7.8 


4.0 


5.1 


(17, 0) 


1.35 


241 


136 


0.028 


6.4 


3.0 


5.0 






904 


SQ 




7 i 


3 7 


3 r > 


(20, 0) 


1.59 


208 


99 


0.027 


6.0 


3.1 


4.5 


(16, 8) 


1.68 


194 


84 


0.026 


6.1 


3.2 


4.1 


(22, 0) 


1.75 


179 


69 


0.025 


6.6 


3.6 


3.3 


(23, 0) 


1.83 


182 


75 


0.026 


5.8 


3.1 


3.9 


(25, 0) 


1.99 


160 


54 


0.024 


6.2 


3.6 


3.8 


(26, 0) 


2.06 


162 


59 


0.025 


5.5 


3.1 


3.7 


(24, 8) 


2.29 


142 


43 


0.025 


5.7 


3.5 


4.2 



TABLE I: Tube indices, diameters in nm, exciton binding energy in meV, and Fo fields in V//im in 
semiconducting carbon nanotubes studied here. The rest of the columns give the best fits for do in nm, 
Ks, Kb, and a in Eq. [2] [3] [4] In italic are shown fit parameters for an averaged fit error larger than 8%. 
Parameters A — 0.009 and j3 = 1.74 are kept fixed in the fits. 

the following equation: 

^and _ . d jedF) 2 
1 XT — -a+ — + k s — -g— [Zj 

1 band i 1 exc w 

where A, do, and k s are fit parameters, d is the tube diameter, F is the field strength, e - electron 
charge, and E\, is the first exciton binding energy. The field dependence in Eq. ([2]) is motivated 
by first order perturbation theory for the bound exciton wavefunction mixed by the field with 
the first band-to-band states separated in energy roughly by the exciton binding energy. The 
matrix element for the coupling is proportional to the field and the exciton size, which in turn is 
proportional to the tube diameter [8(. The best fit for all 13 tubes, studied here, is achieved with 
A = 0.009, do — 0.026 nm, and k s — 6.4. The exact values used for each individual fit for all 13 
tubes are given in table [J In agreement with our predictions the exciton spectral weight reduction 
has been found to scale quadratically with the applied electric field in recent electroabsorption 
measurements [lj| . 
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FIG. 3: Exciton binding energy dependence on electric field in (12, 4) tube with d — 1.15 nm (red circles), 
(16, 0) d = 1.3 nm (blue circles), (19, 0) d = 1.5 nm (green circles), (22, 0) d — 1.7 nm (cyan circles), and 
(26, 0) d — 2.1 nm (magenta circles). The solid curves are the best fits to Eq. © with Kb from table U 



The first exciton binding energy increases quadratically with the field, leading to a red-shift of 
the En absorptionpeak. Despite the fact that the first exciton is nearly degenerate with the lower 
energy dark state [8|, [2CJ, [21 



22| , we do not find enough mixing between these two states to lead 
to an observable linear Stark component. The linear Stark effect would be expected only if the 
off-diagonal coupling strength between the dark and the bright excitons is larger than its energy 
splitting. In zig-zag tubes the mixing is identically zero and it is negligibly small in chiral tubes. 
This is due to the different spatial symmetries of these two states, which lead to a negligible dipole 



moment between them, as was pointed out in Ref [23| ■ The field dependence of the binding energy 
is shown in Fig. [3J It is quadratic at low fields, but it deviates from the parabola at higher fields. 
The quadratic part of the Stark effect can be well described by the following equation: 
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whose form is similar to the second order perturbation theory for the bound exciton energy and 
the exact treatment of the Stark effect in hydrogen atom . A quadratic shift was also observed 
in the quantum confined Stark effect J, 25 1. The best fit for all 13 tubes studied here is achieved 
with Kb = 3.4. Similar arguments apply to the field dependence of the lower energy dark state. As 
a result, we find the splitting between the bright and the dark states to be nearly independent of 
the field strength. 



The second E22 exciton binding energy has a similar quadratic field dependence suggesting a 
red shift of the absorption peak, but the magnitude of the shift is typically smaller by a factor of 
three. The E22 exciton can be coupled by the field to states both higher and lower in energy, which 
would lead to opposite contributions to the sign of the Stark shift. In principle, the sign of the 
net effect may become reversed, depending on the strength of the couplings to the first and second 
band states. For e = 3.3, we find E22 to always red-shift (increase of the binding energy) with 
field, but the net effect is reduced in magnitude. However, when we reduce the dielectric constant 
to e = 2.0, the second exciton binding energy increases, the E22 resonance comes close to the first 
band-to-band absorption peak, and the coupling to the latter dominates leading to an overall E22 
blue shift. At even higher fields of about F > 10 V//im and e = 2.0, we find that the sign of 
the second exciton Stark shift suddenly reverses and we obtain a red-shift. This nonperturbative 
result can be explained by the Franz-Keldysh modulation of the band-to-band states. 
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FIG. 4: Exciton dissociation rate as a function of field in (13, 0) tube with d = 1.0 nm (red), (17, 0) 
d = 1.35 nm (blue), (20, 0) d = 1.6 nm (green), (16, 8) d = 1.7 nm (cyan), and (24, 8) d = 2.3 nm 
(magenta). The black squares are the linewidths from the absorption fit to the Fano lineshape. The solid 
curves are the best fit to Eq. (JTJ with (3 — 1.74 and a from table [I] 



Under an applied field, a bound exciton can dissociate into a free electron and a hole in the first 
continuum band, which may contribute, for example, to photoconductivity 
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dissociation rate determines the bound exciton lifetime and results in the broadening of the En 
absorption peak. At high fields, we can determine the bound exciton dissociation rate by fitting 
the En absorption from the Bethe-Salpeter equation solution to the Fano lineshape. The results of 



the fits are shown in Fig. 0] by the black squares. When the broadening becomes less than 1 meV, 
this procedure can no longer be used. For smaller fields, we determine the lifetime by calculating 
the tunneling probability of the bound exciton into the free electron-hole continuum, from the 
leaking of the exciton wavefunction away from the central region. The exciton wavefunctions for 
different fields are shown in Fig. [5] We calculate the tunneling probability from the wavefunction 
weight at an electron- hole separation larger than 50 nm[27|. To convert the tunneling probability 
into the decay rate, we need to know the "attempt" frequency for tunneling. From the Heisenberg 
uncertainty principle, an averaged momentum of the bound exciton is p — 2nh/\, where A is the 
exciton size. According to the virial theorem, the exciton binding energy Eb is proportional to the 
kinetic energy Ek ~ pv/2, where v is an averaged velocity. In 3D case, Eb — Ek, and, therefore, 
one would expect, in the case of nanotubes, an attempt frequency v/(2X) to scale with Eb/(2i:h). 
By adjusting the attempt frequency to agree with the lifetime calculations from the absorption 
lincwidth, we can obtain the bound exciton dissociation rate in a wide range of electric fields shown 
in Fig. 2J The attempt frequency varies from 1 to 1.5 times Eb/(2irK). The exciton dissociation 
rate can be well-fitted by the following equation: 

r (F)=a£ 6 ^exp(-^ (4) 
where F = fiE^ 1 ' 'm 1 ! 2 / 'eh, m is a reduced exciton mass (m 1 — m e 1 + m h ). This form is 
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motivated by the solution of the hydrogen atom in an applied electric field 
here, that the electron and hole dispersions in CNTs are hyperbolic, unlike parabolic dispersion 
typically used in exciton dissociation description in solids. The best fit for all 13 tubes, studied 
here, is obtained for a = 4.1 and (3 = 1.74. As seen from Fig. 21 at high fields F/Fq > 0.14, the 
tunneling rate from Eq. (TTJ) deviates from the numerical result. The criterion for the full exciton 
dissociation [29| is satisfied at a field strength of F c « Fq/2, It is instructive to compare the present 
results with those obtained for the 2D exciton ionization rates in fields parallel to the quantum 

n 

well plane. [25j. The exciton ionization in the latter is one order of magnitude smaller due to the 
difference in exciton mass and binding energy. 

It has recently been shown that under high bias the energetic carriers in nanotubes excite optical 
phonons 
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32, 



33j | and generate a non-equilibrium phonon distribution 34 , |35j , particularly 
when energy dissipation to the substrate is suppressed as in suspended CNTs [35 1. Therefore, one 
may expect that the produced excitons, for example, in electroluminescence experiments (l^j . can 
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FIG. 5: Optically active exciton wavefunction for different fields F = 5.0 (red), 6.0 (magenta), 7.0 (cyan), 
8.0 (green), and 9.0 (blue) V//im. There is a finite probability for exciton to dissociate, which is found 
from the wavefunction weight outside the central region \r e — r^j > 50 nm. 

experience a phonon assisted dissociation. To explore this mechanism, we need to model the energy 
distribution function of excitons, which is determined by the coupling to the optical phonon bath 
and the spectral function of the pumping source of excitons. The optical phonon bath is assumed 
to be at temperature T op and the pumping source of excitons is assumed to be tuned at the En 
energy with intensity Pq. The bound exciton with zero momentum can absorb an optical phonon 
with momentum — q and be promoted to state q with an optical phonon energy hw p h above the 
bottom of the exciton band E\\. The probability of this event is proportional to the optical phonon 
occupation number n p h(—q) and the strength of the exciton-optical phonon coupling. In the reverse 
process, the high energy g-exciton can emit a phonon and be relaxed to the bottom of the exciton 
band. This probability is proportional to the product of (1 + n p h(—q)) and the strength of the 
exciton-optical phonon coupling. Since the phonon dispersion is much smaller than the exciton 
dispersion, we can reduce the problem to the two level system described by the rate equations: 



dt u v h 



^ = -jv^-^i+M + ^M (5) 

at n T ph T ph 

where Nq and N\ are the ground and vibrationally excited states, n p h is given by the Bose-Einstein 
factor, and the scattering time r p h ~ 30 — 100 fs is determined by the exciton-optical phonon 



coupling strength. The exciton-optical phonon scattering time has been calculated 



and is of 



the same order of magnitude with the measured 
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phonon scattering time. The dissociation rate of the excited state T\ is much faster than that of 
the ground state To. It can be modeled by Eq. Q with binding energy replaced by Eb — fttxiph, 
where hw p h — 200 meV is the optical phonon energy. Note that in the absence of the exciton 
pumping source Po = and the exciton decay channels To^i = 0, equation O gives a steady state 
thermal exciton distribution in equilibrium with the optical phonon sink bath. The net exciton 
dissociation rate, including the phonon assisted mechanism, can be obtained from: 

T(T , _ r iVo + iyvi _ r (i + n ph + T ph r 1 /h) + 

[op) ~ No + m ~ l + 2n ph + T ph T 1 /h U 

where Nq and Ni give the steady state solution in Eq. ([5]). In the limit of a very efficient ther- 
malization, T p hT\ <C 1, the exciton ionization rate from Eq. [6] is given by the thermal average of 
the decay rates of the two states. On the other hand, in the opposite limit T p hTi 3> 1, which is 
realized, for example, when the vibrationally excited state is above the band-to-band continuum, 
Eb < fiwj,;,, the exciton dissociation rate depends also on the exciton-optical phonon scattering 
time as T(T op ) w r + Uph/Tp^. Even at room temperature, the phonon assisted contribution to 
the dissociation rate can be appreciable, especially at low fields, due to the strong exciton-optical 
phonon interaction in carbon nanotubes. In opto-electronic devices operating at high currents, 
the phonon assisted mechanism would dominate the net exciton dissociation rate due to the high 
optical phonon occupation. 

In conclusion, we predict strong modulation of the absorption spectra of carbon nanotubes by 
the electric field. We find a large electric field enhancement of the band-to-band absorption, or the 
Sommerfeld factor, and a large red Stark shift of the bound exciton, especially in large diameter 
tubes. We find that the bound exciton ionization rate depends strongly on field strength and can be 
in the sub-picosecond range at moderately high electric fields. Finally optical phonon occupation, 
which can be high in opto-electronic devices, dramatically increases exciton dissociation. 
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